The pion electromagnetic form factor 
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A ratio of lattice correlation functions is identified from which the pion form factor can be obtained directly. 
Preliminary results from quenched Wilson simulations are presented. 



1. MOTIVATION 

The pion form factor is a convenient way to 
study the transition from perturbative to non- 
perturbative QCD. This transition is typically ex- 
pected to occur at a smaller than for the nu- 
cleon, and should thus be more easily attainable 
in experiments and in lattice studies. Also, the 
asymptotic normalization of the pion form fac- 
tor is known in terms of /,r- In addition, lat- 
tice simulations of the pion form factor bene- 
fit from the complete absence of "disconnected" 
quark diagrams. P 

Some experimental data already exist, another 
experiment is in progress right now at Jefferson 
Lab, and measurements at higher are planned 
for the future. 0. 

Pioneering lattice QCD studies were carried 
out in the 1980's [Tl3l4| . and only recently has 
the pion form factor been revisited using lattice 
methods. |5I6| In light of the substantial exper- 
imental effort, it is important to produce lat- 
tice QCD results with significantly smaller quark 
masses than have yet been employed, and even- 
tually to attain larger values of . Unquenched 
simulations will also be important. Here, after 
a general discussion of correlation functions, the 
issue of smaller quark masses will be addressed 
using quenched Wilson data. 

2. CORRELATION FUNCTIONS 

The electromagnetic form factor is obtained 
in lattice QCD simulations by placing a charged 
pion creation operator at Euclidean time i^, a 
charged pion annihilation operator at i/, and a 
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Figure 1. The quark propagators used to compute 
the pion form factor. 

vector current insertion at t as shown in Fig. ^ 
A standard determination of 12 columns in the in- 
verse of the quark matrix provides the two prop- 
agators that originate from ti. The remaining 
propagator is obtained for any particular annihi- 
lation operator, having definite momentum p/, by 
a two step procedure: the propagator from ti to 

is multiplied by the annihilation operator and 
that product is used as the input for a second 12- 
column linear algebra step. This determines the 
entire thick quark line of Fig.^ The momenta at 
the source (pi) and at the insertion (g) can then 
be varied freely without computing new propaga- 
tors. Since the largest is attained in the Breit 
frame {pf = —Pi), the simple choice of = is 
not always desirable. 

The form factor, F{Q^), is defined by 

(^(p;)|T/,(0)K(K))co„ti„uum (1) 

= Zv {TTipf) \V^iO)\ nip,)) = FiQ^)ip,+pf), 

where V^(x) is the chosen vector current. F{Q^) 
can be extracted from ratios of lattice correlation 
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functions. The three-point correlator is 
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where A £ {L,S) and B e {L,S) denote either 
"local" or "smeared". Inserting complete sets of 
hadron states and requiring <tit<^tf, gives 
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Similarly for the two-point correlator, 

T^^iU,tf,p)~.{0\(t>B{x^)\7Tip)) 

x(4p)|^t^(x,)|o)^e-(*/-*')^. (4) 

Finally, consider operators that are temporally lo- 
cal, but spatially smeared democratically among 
the three spatial axes: 



a^{0\Mx)\HP)) = Zs{\p\)e"-P. 
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From this discussion, it is found that the following 
expression for the form factor is independent of 
Zl, Zs{\p\), and all time exponentials: 



T^^{U,t,p,)TCB(^t,^tf,pf) 

I 2ZyEf 
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where the indices A, B and C can be either L 
(local) or S (smeared). For a conserved current, 
Zv = 1. 

3. ACTION AND PARAMETERS 

Our initial explorations have used the 
(quenched) Wilson action with /3 = 6, so a ~ 0.10 
fm. Periodic boundary conditions are used ex- 
cept for the quark temporal boundaries, which 
are Dirichlet. Hopping parameters and lattice 
sizes are shown in Table ^ the resulting pseu- 
doscalar and vector meson masses are also listed. 



Table 1 

The chosen hopping parameters and numbers 
of lattice sites, along with the computed pseu- 
doscalar and vector meson masses. 
K lattice amps amy 



0.1480 16^ X 32 

0.1520 16^ X 32 

0.1540 16^ X 32 

0.1555 24^ X 32 

0.1563 243 X 32 

0.1566 243 X 32 & 32^ x 48 



0.673 0.712 

0.477 0.549 

0.364 0.468 

0.259 0.398 

0.179 0.358 

0.145 0.343 



Notice that the smallest pion mass is less than 
300 MeV, and therefore significantly below those 
used in previous studies of the pion form fac- 
tor. The pion operators are fixed at = 7 and 
t/ = 21 except for the 32"^ x 48 lattices, where 
tf ~ 28. The standard point-split conserved vec- 
tor current is used. A gauge-covariant Gaussian 
smearing is employed at the source. 
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1 + ^ h{x). 
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The various results shown here are based on the 
analysis of between 80 and 155 configurations. 

4. PRELIMINARY RESULTS 

Representative plots of the data are shown in 
Fig. 121 For each of the smaller hopping param- 
eters, there is a reasonable plateau at timesteps 
sufficiently far from the pion operators (timesteps 
and 14). For k = 0.1563 and 0.1566, the possi- 
ble plateau spans only a few timesteps. 

Correlated fits to the data have produced Fig.|3| 
which shows the form factor as a function of 
for each available k. For comparison, the 
monopole curve 1/(1 + Q'^/m^) is also shown 
where nip denotes the physical p meson mass. 
Notice that the experimental measurements from 
Ref. 7 lie on the monopole curve. For each k 
the lattice results also display a monopole form, 
and F{Q^) is found to decrease with decreasing 
quark mass as expected since the effective vector 
meson mass of the monopole form also shrinks 
with decreasing quark mass. 

In conclusion, Eq. {Tj) is found to be a useful 
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Figure 2. Pion form factor data versus timestep, for (a) k=0.1520, (b) k=0.1540 and (c) k=0.1563. 
Numerical values of in GeV are shown in the legends. 
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Figure 3. Results for the pion form factor as a function of for each of the available k values. Experi- 
mental measurements 7 and the monopole approximation are also shown. 



method for extracting the pion form factor. Sim- 
ulations with domain- wall fermions are underway. 
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